Higher-order anisotropic flows and dihadron correlations in Pb-Pb collisions at 
^snn = 2.76 TeV in a multiphase transport model 
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Using a multiphase transport model that includes both initial partonic and final hadronic scatter- 
ings, we have studied higher-order anisotropic flows as well as dihadron correlations as functions of 
pseudorapidity and azimuthal angular differences between trigger and associated particles in Pb-Pb 
collisions at ^snn = 2.76 TeV. With parameters in the model determined previously from fitting 
the measured multiplicity density of mid-pseudorapidity charged particles in central collisions and 
their elliptic flow in mid-central collisions, the calculated higher-order anisotropic flows from the 
two-particle cumulant method reproduce approximately those measured by the ALICE Collabora- 
tion, except at small centralities where they are slightly overestimated. Similar to experimental 
results, the two-dimensional dihadron correlations at most central collisions show a ridge structure 
at the near side and a broad structure at the away side. The short- and long-range dihadron az- 
imuthal correlations, corresponding to small and large pseudorapidity differences, respectively, are 
studied for triggering particles with different transverse momenta and are found to be qualitatively 
consistent with experimental results from the CMS Collaboration. The relation between the short- 
range and long-range dihadron correlations with that induced by back-to-back jet pairs produced 
from initial hard collisions is also discussed. 

PACS numbers: 25.75.-q, 12.38.Mh, 24.10.Lx 
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I. INTRODUCTION 



Heavy ion collisions at the Rclativistic Heavy-Ion Col- 
lider (RHIC) have made it possible to obtain important 
information on the properties of the quark-gluon plasma 
(QGP) that is produced in these collisions [lHll- Among 
the many observables that have been used to infer the 
properties of the QGP are the anisotropic flows of pro- 
duced particles, particularly the elliptic flow (i> 2 ) [H that 
corresponds to the second-order harmonic coefficient of 
the anisotropic flow and is appreciable in noncentral col- 
lisions of identical nuclei. More recently, it was found 
that the triangular flow 113 that corresponds to the third- 
order harmonic coefficient of the anisotropic flow could 
also be significant in collisions of all centralities as a re- 
sult of initial density fluctuations in the collision geom- 
etry 0, 0). Furthermore, it was shown that studying V3 
together with V2 can help better constrain the proper- 
ties of the QGP [1, Also, it was suggested that use- 
ful information about the QGP could be obtained from 
the dihadron correlation [lfj between lower-py associ- 
ated particles from the medium response to the away- 
side jet and higher-py trigger particles produced from 
the near-side jet. Studies based on the transport model 
showed that higher-order anisotropic flows, especially the 
triangular flow, had a large effect on dihadron correla- 
tions induced by initial energetic jets [HI, [Til- To extract 



from the dihadron correlation information on the inter- 
actions of energetic jets produced from initial hard col- 
lisions with medium partons in the quark-gluon plasma 
requires, however, the subtraction of the background con- 
tribution from anisotropic flows and their fluctuations, 
and this is a topic of great current interest 0, [lll - [l^ | . 

Recently, results from the first Pb-Pb collisions at 
t/sJ^ = 2.76 TeV at the Large Hadron_Collider (LHC) 



have attracted a lot of attentions |15H18| . In the most 
central collisions (0 — 5%), the multiplicity density of 
produced charged particles at mid-pseudorapidity is 2.2 
times of that in Au+Au collisions at ^s NN = 200 GeV 
at RHIC. The measured elliptic flow, which can provide 
more reliable information about the produced QGP [lj| , 
is of a similar magnitude as that measured at RHIC. Be- 
sides, experimental data on higher-order flows EU and 
dihadron correlations 
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have also become available, 
providing the opportunity to study the properties of the 
hotter and denser QGP formed at LHC. In particular, it 
was found that the long-range dihadron azimuthal corre- 
lations could be entirely accounted for by the anisotropic 

flows mmi. 

In a recent work [22| . we studied Pb-Pb collisions at 
y/s NN = 2.76 TeV in a multiphase transport (AMPT) 
model. The AMPT model is a hybrid model with the 
initial particle distribution generated by the heavy ion 
jet interaction generator (HIJING) model [23j . In the 
version of string melting, which was used in the pre- 
vious work [22j and will be used in the present study, 
all hadrons produced in the HIJING model through 
Lund string fragmentation are converted to their valence 
quarks and antiquarks, whose evolution in time and space 
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is modeled by Zhang's parton cascade (ZPC) model [241 ] . 
After their scatterings, quarks and aniquarks are con- 
verted via a spatial coalescence model to hadrons, and 
the scatterings among them are described by a relativis- 
tic transport (ART) model [25j]. For a recent review of 
the AMPT model, we refer the readers to Ref. [H] . By 
adjusting the values of parameters in the Lund string 
fragmentation function and the parton scattering cross 
section, we fitted the differential elliptic flow at 40 — 50% 
centrality window and reproduced reasonably well the 
centrality dependence of both the multiplicity density 
and the elliptic flow of mid-pscudorapidity charged par- 
ticles }22| , although the transverse momentum spectrum 
at most central collisions is softer than the data from the 
ALICE Collaboration as a result of the use of massless 
partons in the model that are less affected by the radial 
flow [261 ] . In the present paper, we study higher-order 
anisotropic flows and dihadron correlations in Pb-Pb col- 
lisions at LHC by using the AMPT model with the newly 
fitted parameters. Our results show that in heavy ion col- 
lisions at LHC higher-order flows are appreciable and a 
ridge structure exists in the dihadron correlation as seen 
in the experiments |2l| . The ridge structure is a result 



of the longitudinal expansion of the anisotropic flows as 
pointed out previously in Ref. [13 ■ 

This paper is organized as follows. We first give some 
general discussions on anisotropic flows and dihadron cor- 
relations in Sec. |TIJ We then show in Sec. IIII Al the re- 
sults on the centrality dependence of anisotropic flows 
and their transverse momentum dependence in most cen- 
tral and mid-central collisions and in Sec. IIII Bl those on 
the dihadron correlations in most central collisions. A 
summary is then given in Sec. IIV1 



II. ANISOTROPIC FLOWS AND DIHADRON 
CORRELATIONS 

To help understand the results of anisotropic flows and 
dihadron correlations in the following sections, we first 
discuss generally the contribution from anisotropic flows 
and their fluctuations as well as the nonflow contribution 
to the two-dimensional dihadron correlation between two 
particles. 

The momentum distribution of produced particles in a 
heavy ion collision event can be generally written as 



the orthonormal relation of the harmonic terms, i.e. 
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By neglecting the correlation between particle number 
and flow, the event average of the two-particle correlation 
can then be approximated by 
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where (• ■ ■ ) e denotes average over all events. 

Assuming that the dependence of the anisotropic flows 
on the pseudorapidity is weak, the second term in Eq. ([3]) 
can be taken out of the integration and further expressed 
approximately as 

(«n(Pr, il) v n(PT, V + A?/)) cos(nA</>) 

~ (v n (PT))c(Vn{PT))cCOs(nA<f>) 

+ FF{vM),v n (p b T )) cos(nA^) + NF(A& An),(4) 

where FF(v n (p^,),v n (p b r )) is the flow fluctuation contri- 
bution and NF(A<f>, An) is the non-flow contribution due 
to jet correlations or resonance decays that are impor- 
tant only for small An and A<p ~ or x. For the case of 
Pt = p\i t ne usual anisotropic flows calculated using the 
two-particle cumulant method [1^, [2t| after integrating 
the pseudorapidity difference is then obtained, i.e., 



v n {2} = V(cos(nA0)). 
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where <j) is the azimuthal angle, ^„ is the nth-order event 
plane angle, and N(pt,t)) and v n (pT,r]) are the num- 
ber of particles of transverse momentum px and pseu- 
dorapidity rj and their nth-order anisotropic flows, re- 
spectively. The correlation between two particles with 
transverse momenta p?p and p^ and azimuthal angular 
difference A<j) and pseudorapidity difference An within 
the same event can be calculated from Eq. (JT|) by using 



Using the fact that the first term in Eq. ([3]) can be 
calculated from the mix-event correlation, i.e., 
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since the harmonic terms vanish due to the independent 
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event plane angles in different events, it follows that 

/ d 2 N£™ \ / cPNg? \ 
\ dArjdAcj) ] 1 \ dAr]dA(f> J 

« 1 + (w„(pT))e(fn(PT))eCOs(nA0) 

+fFK(PT)> «n(Pr)) cos(nA0) + AF(A0, At?). 

(7) 

The dihadron correlations thus have contributions from 
anisotropic flows and their fluctuations as well as the non- 
flow contribution. 



III. RESULTS 
A. Anisotropic flows 



they are slightly larger at small centralitics. This is due 
to the fact that we have chosen in our previous work (22[ 
to fit the elliptic flow at the 40 — 50% centrality window. 
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FIG. 1: (Color online) Centrality dependence of v n (n = 
2,3,4) for mid-pseudorapidity (\r]\ < 0.8) charged particles 
obtained from the two-particle cumulant method in Pb-Pb 
collisions at y/s NN = 2.76 TeV from the string melting AMPT 
model. The ALICE data (filled symbols) are taken from 
Ref. [ITI. 

We first consider the centrality dependence of 
anisotropic flows v„(n = 2,3,4) for mid-pseudorapidity 
charged particles in Pb-Pb collisions at ^/s NN = 2.76 
TeV from the string melting AMPT model using the 
two-particle cumulant method (Eq. ([5])). To reduce non- 
flow effects on the calculated v n , only pairs with pseu- 
dorapidity difference \Arj\ > 1 are considered in the 
present study as in experimental analyses pol ] . Results 
for transverse-momentum-integrated (0.2 < px < 5.0 
GeV/c) anisotropic flows are shown by open symbols 
in Fig. [Tl Compared with the ALICE data (filled sym- 
bols) [20] , the AMPT model reproduces reasonably well 
the centrality dependence of v n (n = 2,3,4) except that 



FIG. 2: (Color online) Transverse momentum dependence 
of vi (a), «3 (b), V4 (c), and vs (d) for mid-pseudorapidity 
(\r]\ < 0.8) charged particles obtained from the two-particle 
cumulant method in most central (0 — 5%) Pb-Pb collisions 
at y/s^ = 2.76 TeV from the string melting AMPT model. 
The ALICE data (squares) are taken from Ref. [2(|]. Results 
for | At? | > 0.2 and |A^| > I are shown by filled and open 
symbols, respectively. 
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FIG. 3: (Color online) Same as Fig. [2] but for centrality 30 — 
40%. 

In Figs. [2] and O we display by circles the trans- 
verse momentum dependence of anisotropic flows v„ (n = 
2,3,4,5) of mid-pseudorapidity (|?7| < 0.8) charged par- 
ticles in the most central (0 — 5%) and the mid-central 
(30 - 40%) Pb-Pb collisions at = 2.76 TeV, re- 
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spectively, from the AMPT model with string melting, 
and compare them with the ALICE data (squares) from 
Ref. ■ It is seen that the anisotropic flows for different 
pscudorapidity cut \Ar}\ > 0.2 and IA77I > 1 are similar 
to those from the ALICE data, while at larger transverse 
momenta in the case of the smaller pseudorapdity cut the 
AMPT model gives larger values due to overestimated 
nonfiow effects. As for RHIC [l2j, anisotropic flows for 
the most central collisions from the AMPT model are 
overestimated in comparison with those from the ALICE 
data, while they are consistent with the ALICE data for 
mid-central collisions. The reason for this might be due 
to the constant parton scattering cross section used in 
the AMPT model for all centralities. Since the temper- 
ature of produced partonic matter is higher in central 
collisions than in mid-central collisions, a smaller str ong 
coupling constant and a larger screening mass [3(| Ef 
should be used in calculating the parton scattering cross 
section. This is expected to lead to a smaller parton 
scattering cross section and smaller anisotropic flows for 
central collisions. We note that once the elliptic flow is 
fitted, the higher-order flows are automatically consis- 
tent with the experimental data. This indicates that the 
HIJING model, as an initial distribution generator, pro- 
vides the correct initial spatial anisotropies through the 
Glauber model calculation. 




FIG. 4: (Color online) Two-dimensional dihadron correla- 
tion per trigger particle as a function of Arj and A<f> for 
1 < pf BOC < 2 GeV/c and 2 < < 3 GeV/c for - 5% 
most central Pb-Pb collisions at ^/s NN — 2.76 TeV from the 
string melting AMPT model. 



Dihadron correlations 



As in Ref. [2l| 
according to 



we calculate the dihadron correlation 
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is the raw correlation per trigger particle from pairs in 
same events, 
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is the background correlation per trigger particle from 
pairs in different events, and B(Ar] = 0, A0 = 0) is the 
normalization factor. Eq. ([SJ is thus similar to Eq. (J7J in 
Sec. in 

In Fig. 2] we show the two-dimensional dihadron cor- 
relation per trigger particle from the string melting 
AMPT model in the - 5% most central Pb-Pb colli- 



sions 



at ^ft 



NN 



2.76 TeV for trigger particles and as- 



sociated particles in the transverse momentum windows 
2 < pr is < 3 GeV/c and 1 < p^ soc < 2 GeV/c, respec- 
tively Similar to the experimental results in Ref. |21| . 
there is a peak at A77 = and a ridge structure extend- 
ing to I At; I = 4 at the near side (A</> ~ 0) as well as 
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FIG. 5: (Color online) Same as Fig. [4] but without final-state 
interactions (FSI). 



a broad structure at the away side which also extends 
to | A 77 1 = 4. The ridge structure was first discovered in 
Ref . I32l . Many explanations have since been proposed 
[33l - l39l ] until it was recently realized that the ridge was 
dominated by the anisotropic flows jlil . I20I [2lj as can 
be seen from Eq. ([?])• To illustrate this effect, we re- 
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peat the calculation by turning off both partonic and 
hadronic scatterings in the AMPT model, and the re- 
sults are shown in Fig. [5] It is seen that the near-side 
ridge now disappears and the away-side broad structure 
is also largely weakened as both are results of the col- 
lective flow that is generated by final-state interactions 
(FSI). Furthermore, the peak along the Arj direction is 
sharper than that in the case with FSI. We note that the 
wider nonfiow contribution in Arj direction in the case 
with FSI is related to the early explanations of the ridge 
structure HHH. 




by taking average of the two-dimensional dihadron corre- 
lation over the corresponding Arj windows plj as shown 
in the left and middle columns of Fig. [7] for different trans- 
verse momenta of trigger particles. Compared with the 
short-range correlations, the long-range correlations have 
much weaker peaks at the near side due to the absence 
of the nonfiow contribution. The strength of the near- 
side peaks, i.e., the near-side associated yield, decreases 
with increasing px of trigger particles for both short- and 
long-range correlations, although their heights are simi- 
lar for short-range correlations. In addition, the strength 
of away-side correlations decreases with increasing px of 
trigger particles for both short-range and long-range cor- 
relations. 
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FIG. 6: (Color online) The near-side associated yield as a 
function of |Ar;| (a) and those as functions of trigger parti- 
cle transverse momentum in the jet region (b) and the ridge 
region (c) for — 5% most central Pb-Pb collisions from the 
string melting AMPT model. 

Above results can be quantified by displaying the \Arj\ 
dependence of the near-side associated yield. Taking the 
near side as \A<p\ < A</>zyam, with A(/>zyam obtained 
using the zero-yield-at-minimum (ZYAM) method poj to 
fit the dihadron azimuthal correlation by a second-order 
polynomial in the range 0.5 < A<f) < 1.5 as in Ref. plj . 
the results are shown in Fig. [5] (a) for different transverse 
momentum windows of trigger particles (2 < p^. Ig < 3 
GeV/c, 3 < pr is < 4 GeV/c, and 4 < P ^ is < 5 GeV/c) 
but the same transverse momentum window of associated 
particles (1 < p^ ssoc < 2 Gev/c). It shows a similar 
decreasing trend with increasing |A?y| as in Ref. [2l|, and 
it is broader for the case with lower-p^ trigger particles. 
The width of the near-side associated yield is related to 
the broadening of the jet cone affected by the longitudinal 
flow, and it seems that the jet cone from higher-pr jets 
is more focus and less affected by the longitudinal flow. 

For the dependence of the near-side associated yield 
on the trigger particle transverse momentum, we first 
introduce the short-range (|Arj| < 1) and long-range 
(2 < \Ar]\ < 4) dihadron azimuthal angular correlations 



FIG. 7: (color online) Short-range (\Ar)\ < 1, (a), (d) and 
(g)) and long-range (2 < |A7j| < 4, (b), (e) and (h)) dihadron 
correlations per trigger particle as functions of the azimuthal 
angular difference A<f> for different transverse momentum win- 
dows of trigger particles from — 5% most central Pb-Pb colli- 
sions in the string melting AMPT model with (solid lines) and 
without (dashed lines) final-state interactions. The difference 
between the short-range and long-range dihadron azimuthal 
correlation is shown in panels (c), (f), and (i). 

Results on the trigger particle transverse momentum 
dependence of the near-side associated yield in the jet 
region and the ridge region are shown in Fig. [5] (b) and 
(c), respectively. In the ridge region, i.e., the long-range 
correlation (2 < \Arj\ < 4), the near-side associated yield 
decreases with increasing px of trigger particles. This is 
understandable as the near-side ridge is dominated by the 
collective flow that does not affect much higher-p^ par- 
ticles. The near-side associated yield in the jet region, 
which is the difference between the associated yield in 
the long-range correlation (2 < \Arf\ < 4) and the short- 
range correlation (|A?y| < 1), increases with increasing 
Pt of trigger particles. This can be again understood as 
higher-j»T particles are more likely to be accompanied by 
a larger number of lower-p^ particles because they can 
fragment into or dump energy to more lower-py particles. 
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All these a re q ualitatively consistent with the experimen- 
tal results [Hf. 

Since the dihadron correlation includes contributions 
from both the anisotropic flows and the medium response 
to the back-to-back jet pairs produced in initial hard col- 
lisions, to isolate the latter effect requires the subtraction 
of the flow contribution. As discussed in Refs. fl4l . l2Cil . l2lj ]. 
the long-range dihadron correlation is largely caused by 
the anisotropic flows and their fluctuations. We can 
therefore use the long-range dihadron azimuthal corre- 
lation shown in the middle column of Fig^ 7] as an esti- 
mation of the background contribution [l4| by neglect- 
ing the pseudorapidity dependence of the anisotropic 
flows and possible non-flow contributions for larger |A?7|. 
The resulting background-subtracted correlations, i.e., 
the difference between the short-range correlations and 
the long-range correlations, are shown in Fig. [7] (c), (f), 
and (i). It is seen that the resulting away-side correla- 
tions are very weak for all cases and this is similar to 
the results in Ref. [IH obtained using a different method 
to estimate the background correlations for central colli- 
sions. In addition, the near-side correlations in this case 
are stronger for highcr-p^ trigger particles, which can be 
understood as the stronger medium response to higher-pr 
jets. For comparisons, we also show in Fig. [7] by dashed 
lines the dihadron azimuthal correlations obtained with- 
out final-state interactions, which are similar to those 
from p+p collisions at the same energy. The background- 
subtracted away-side correlations in this case is similar to 
that obtained with final-state interactions. The smear- 
ing of the away-side jet in heavy ion collisions without 
final-state interactions or p+p collisions is due to the 
fairly low px of trigger and associated particles consid- 
ered here, which can lead to an away-side jet with very 
different pseudorapidity from that of the triggered jet [41j 
as shown in the middle column of Fig. [Jj A peak in the 
away-side correlations of such collisions would, however, 
appear if the associated particles have similar momenta 
as those of the trigger particles or the trigger particles 
have very high pt- This is very different from the case 
where final-state interactions are included. Because of jet 



quenching as a result of final-state interactions, no peak 
structure is seen in the away-side correlations unless the 
trigger particles have extremely high pt- 



IV. CONCLUSIONS AND DISCUSSIONS 

We have studied higher-order anisotropic flows and di- 
hadron correlations in Pb-Pb collisions at ^/s NN = 2.76 
TeV within a multiphase transport model with parame- 
ters fitted to reproduce the measured multiplicity density 
of mid-pseudorapidity charged particles in central colli- 
sions and their elliptic flow in mid-central collisions in 
the previous work. We have found that the resulting 
higher-order anisotropic flows slightly overestimate the 
experimental data at small centralities but are consis- 
tent with them at other centralities. We have obtained 
the ridge structure along the pseudorapidity direction in 
the near side of the dihadron correlations, and it disap- 
pears when final-state interactions are turned off in our 
model. We have also studied both the short-range and 
long-range dihadron azimuthal correlations for different 
transverse momenta of trigger particles, and they are seen 
to be quantitatively consistent with experimental results. 
We have further attempted to determine the background- 
subtracted short-range dihadron azimuthal correlations 
by taking the long-range dihadron azimuthal correlations 
as the background, and they are found to be similar to 
those obtained previously using a different method. 
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